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ABSTRACT 

One of the challenges that models of AGN heating of the intracluster medium (ICM) face, is 
the question how the mechanical luminosity of the AGN is tuned to the radiative losses of the 
ICM. Here we implement a simple ID model of a feedback mechanism that links the luminosity 
of the AGN to the accretion rate. We demonstrate how this simple feedback mechanism leads 
to a quasi-steady state for a broad range of parameters. Moreover, within this feedback model, 
we investigate the effect of thermal conduction and find that its relative importance depends 
strongly on the cluster mass. 

Subject headings: galaxies: clusters: general, galaxies: active, X-ray: galaxies: clusters 



1. Introduction 

Currently, the most popular model that is in- 
voked to explain the dearth of gas below about 1 
keV in the ICM relies on the heating by a cen- 
tral AGN (Binney & Tabor 1995, Tabor & Binney 
1993, Churazov et al. 2001, Briiggen & Kaiser 2002 
, Briiggen et al. 2002 Reynolds et al. 2001). Nu- 
merous observations of x-ray deficient bubbles in 
clusters and of motions induced by these bubbles 
have substantiated this model (e.g. Fabian et al. 
2003, Mazzotta et al. 2002, Saxton et al. 2001, Mc- 
Namara et al. 2001, Blanton et al. 2001). On the 
theoretical side, the efficiencies with which bub- 
bles heat the ICM has been investigated, both nu- 
merically and analytically. E.g. Ruszkowski et al. 
(2004) found that up to 50 % of the internal en- 
ergy of these bubbles is dissipated through viscous 
dissipation of waves and pdV work. The heating 
profile caused by bubbles has been parametrized 
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in a prescription named 'effervescent heating' by 
Begelman (2001), which will be used later in this 
paper. 

However, one of the challenges that models of 
AGN heating face, is the question how the me- 
chanical luminosity of the AGN is fine-tuned to 
the radiative losses of the ICM. If AGN heating is 
responsible for a diminution of the mass deposition 
rate in a whole range of clusters, the AGN needs 
to be regulated by the ICM itself. It has been sug- 
gested by some authors that AGN feedback may 
play a crucial role in self-regulating cooling flows 
(e.g. Churazov et al. 2002, Ruszkowski & Begel- 
man 2002, Brighenti & Mathews 2003). 

In this paper we present a simple model by 
which the central AGN adjusts its luminosity. 
This model assumes that the mechanical luminos- 
ity of the AGN is proportional to the accretion 
rate onto the AGN. In this picture, a very lumi- 
nous central source will lead to a decrease in the 
accretion rate and hence also to a decline of the 
luminosity until, eventually, the cluster reaches a 
quasi-steady state. With the aid of a spherically 
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symmetric ID model wc study under which cir- 
cumstances a self-regulated cooling flow is estab- 
lished. For feedback to be a viable model, the pa- 
rameters of this model should not need any tuning. 

In the following we will describe our initial mod- 
els, our computational method and the implemen- 
tation of feedback in more detail. Then, in section 
4, we will present and discuss our results. 

2. Model 

2.1. Initial profiles 

The initial cluster profiles are computed within 
the present canonical cosmological parameters, see 
Tab. 1. For the dark matter density distribution 
we assume an NFW profile given by (Navarro et al. 
1996) 

Pdm = ScT''^ {rs +r)~^, (1) 

where 6c is the characteristical density parame- 
ter and Tg the scaling radius. Thus, a cluster 
model is characterized by the virial mass Mvir 
and by the concentration parameter c = Ts/r^ir- 
Note, that the virial radius can be derived from 
4/3 Trr^i^ h? p„it^c = M^-„, where Pcrit is the crit- 
ical density and the density contrast Ac can be 
approximated by Ac = 178 ng '*^ (Eke et al. 1998). 

Unlike for dark matter distributions, there is 
no universal profile for the intracluster gas. Here, 
we follow Roychowdhury et al. (2004) and assume 
that the initial temperature distribution is given 
by (Loken et al. (2002)) 

(2) 



TlnitiaK?") = 1.3 Tew (1 + 1.5r/rvir) 



Now, the emission- weighted X-ray temperature. 
Tew, is known to scale with the cluster mass. Here, 

we adopt the Tx — M-relation by Sanderson et al. 
(2003) who have analysed a large sample of virial- 
ized systems: 

M200 = 2.34 X IO^^Mq (j^) . (3) 

Moreover, we assume that the gas is in hydrostatic 
equilibrium, i.e. the pressure force is balanced by 
gravitation 

1 dP GMtotai(<r) 



Pgas dr 



(4) 



where G is the gravitational constant and Mtotai(< 
r) is the total cumulative mass up to r. We can 
express the pressure as P = nksT, where n is 
the particle number density, which is given for 
a fully ionized gas of primordial composition by 

n = nnil + riHelll + «e = nH(8-5lHe)/(4-4lHe), 

where Yue is the Helium mass fraction. Since 
we assume an initial temperature profile, we can 
rewrite Eq. (4) and solve for the Hydrogen density 
riH 



dnu n-R ( nmn 1 x , dT 



(5) 



where the mean molecular weight is given by /x = 
4/(8 - SFhc)- We integrate Eq. (5) subject to the 
condition that at radius r2oo the ratio of the cu- 
mulative masses Mg^s/Mdm is equal to the cosmo- 
logical value riB/i^Q — ris). 



2.2. Time integration 

The entropy index a = T/ n'J'^ at a given radius 
evolves with time according to 



dfT(r) 
At 



^ ^ 2 1 
3 P(r) 



Mr)-r(r)}, (6) 



where W(r) denotes the total heating and r(r) the 
cooling cooling rate. For the subsequent calcula- 
tions it is convenient to express the entropy pro- 
file as a function of the total enclosed mass, i.e. 
a = a{M). Starting from an initial entropy pro- 
file that is calculated from the models described 
above, the cluster is evolved in time. 

After each time step. At, the new profiles are 
calculated in two steps: first we calculate the 
change of entropy subject to heating and cooling, 
i.e. 

da{r) 



dt 



-At . 



(7) 



In the second step we calculate the new hydro- 
static equilibrium into which the halo settles by 
solving 



dr 
dM 
dP 
dM 



1 



1 



47rr^Pgas(?') 47rr^ 

gMtotal(< r) 

47rr^ 



f) 



1/7 



(8) 
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Equations (8) are discretized on a non-uniform 
grid. A fine grid covers the inner region of the 
cluster out to a mass that in the initial model cor- 
responds to a radius of 5 x 10~^rvir and a coarser 
grid that covers the outer region out to r2oo- Thus, 
the mass of a shell in the finely resolved region is 
0.1% of the cumulative mass Mga,s(< r) compared 
to 2% in the coarse region. This results in about 
4000 mass shells altogether. 

Equation (8) is solved subject to the boundary 
condition that the pressure at r2oo remains con- 
stant during the entire evolution, i.e. P{r2oo) = 

-P200- In practice, we solve the system of coupled 
differential equations (8) subject to the initial con- 
ditions: 



M{0) = 



(9) 



where Pq is a trial value which is varied until the 
pressure at r2oo is equal to the specified pressure 
P200- Note that this implies that the radius r2oo 
changes in the course of time as the cluster ex- 
pands and contracts. 

2.3. Radiative cooling and conduction 

In the temperature range that we are mainly in- 
terested in, namely between 10^ — 10* K, radiative 
cooling is dominated by free-free emission. For a 
fully ionized plasma the cooling rate can be ap- 
proximated by (Katz et al. 1996) 



Afree-free = 1.42x 10-2VT^/^(nH+4nHe)ne 



erg 

s cm^ 
(10) 



with the Gaunt factor 

gs = 1.1 + 0.34exp{-(5.5 - logT)V3}. (11) 

In addition we consider metal-dcpcndcnt cooling, 
which we approximate crudely following Theis 
et al. (1992) 



A. 



metal 



2Q-22-5m(X)-F7VZ ym(Z) ^2^ "^^g 



scm'^ 



m{Z) = 



2.5 + 7y/Z 



:,(12) 



5 - log(1.48 X 1011^11 + 106) ' 

where Z denotes the mass fraction of metals. Here 
we assume Z = 0.3 Zq and we limit the cooling 



in the core by assuming a minimum entropy index 

Cmin = 10^^ kcVcm^. 

The role of thermal conduction in the ICM has 
been the subject of a long debate and, owing to the 
complex physics of MHD turbulence, the value of 
the effective conductivity remains uncertain. The 
thermal conductivity of an unmagnetised, fully 
ionised plasma was calculated by Spitzer (1962). 
Originally it was thought that the magnetic field 
in clusters strongly supresscd the thermal conduc- 
tivity because the magnetic fields prevented an ef- 
ficient transport perpendicular to the field lines. 
Even if the transport can be efficient along the 
magnetic field lines, the overall isotropic conduc- 
tivity was thought to be many orders of magnitude 
less than the Spitzer value. This paradigm has 
been supported by a number of observations, such 
as sharp edges at so-called cold fronts and small- 
scale temperature variations (Markevitch et al. 
2000, Vikhlinin et al. 2001). 

Recent theoretical work by Narayan & Medvedev 
(2001), Chandran et al. (1999), Chandran & Cow- 
ley (1998) and earlier work by Rechester & Rosen- 
bluth (1978) have shown that a turbulent mag- 
netic field is not as efficient in suppressing thermal 
conduction as previously thought. It is argued 
that chaotic transverse motions of the tangled 
magnetic field lines can enhance the cross-field 
diffusion to an extent that the effective conduc- 
tivity is of the order of the Spitzer value. Thus, 
we decided to investigate the effects of thermal 
conduction within our feedback model. Recently, 
the role of thermal conduction in clusters was also 
investigated numerically by Dolag et al. (2004) 
and Jubelgas et al. (2004). 

The energy fiux due to thermal conduction is 
given by 

i^cond = -KVr, (13) 

where k is the coefficient of thermal conductivity 
which we assume to be a fraction / of the Spitzer 
conductivity 

KoT^/^ ~ 5 X 10-^ . (14) 

scmK 

Under the assumption of spherical symmetry, the 
heating and cooling rate due to thermal conduc- 
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tion become 



■ cond 



^2 (^j. 



{r'f'cond} 



-/spitzKo <! 2T— + -r 



dr 



dr 



The timestepping in the integration of Eq. (8) 
is constrained by the appropriate Courant condi- 
tion (Ruszkowski & Begelman 2002). The most 
restrictive term comes from thermal conduction 
which can be obtained from a von Neumann sta- 
bihty analysis, i.e. 

M^^h^^^, (16) 

where Ar denotes the local spatial resolution. The 
time step obviously changes in the course of the 
simulation but for our resolution a typical time 
step is of the order of 0.05 Myr. 

Since we do not use a uniform grid, especially 
with respect to the spatial spacing, we compute 
the derivatives of the temperature by a local poly- 
nomial fit. Around each point in radius we expand 
the temperature to second order 

HT 1 d^T 

T{r) = T{ro) + ^(r- ro) + ^^(r-rof (17) 

and approximate the derivatives with a weighted 
least-square regression. We solve the equation 

/ T{ro) \ 
Q = W • dT/dr (18) 
V l/2d2T/dr2 / 

where 

Wij = S„a;„(r„-ro)'+^'. (19) 

The sum is taken over all sampling points n and 
i,j e {0, 1, 2}. The weights are given by 

= I {rn - r„-i) e-4('-"-'^'')VA^ : |r„ - ro\ < 
[ : otherwise 

(20) 

For the smoothing length we take A = max{0.05 x 
r, O.OSkpc}. Finally, we neglected convection 
which proved to be virtually irrelevant for our sim- 
ulations since the entropy remained monotonous 
even in the presence of strong heating. 



2.4. Heating and feedback 

(15jn the centres of clusters active galactic nuclei 

inflate buoyant bubbles of relativistic gas which 
cPT 1 

_l_ rT- rs^ease some of their internal energy by pdV -work 

^''^oji the ambient medium. This work is done as 
the bubbles expand on their ascent through the 
stratified cluster medium. Churazov et al. (2001) 
and Begelman (2001) have calculated this work 
and found that the heat released by one bubble is 
given by 



h{r) oc p(T''-i)/T' 



1 d InP 



r d Inr 



(21) 



Under the assumption of spherical symmetry the 
heat has to be distributed over spherical shells, 
which introduces a geometric factor /geom oc 
1 / 47rr^. However, at radii smaller than a 
typical bubble diameter this argument does not 
hold. In the center we assume that the heating 
is given by the heating function of one bubble, 
therefore we cut off the geometrical factor below 
the typical bubble size ^bubble, by changing /geom: 
— > + ^bubble- -^^ bubbles rise through the 
ICM they are likely to be disrupted by Rayleigh- 
Taylor instabilities. The smaller fragments rise 
more slowly and are disrupted more easily until 
they virtually stop. As a result the effervescent 
heating rate is effectively cut off at the periph- 
ery of the cluster. This effect is approximated by 
introducing a factor exp{— (r/rdisr)}) where the 
disruption radius, r^isr, is set to 0.5 Mpc. We 
note, however, that this cut-off at large radii only 
has a mild effect on our simulations and is not 
strictly necessary. Thus the effervescent heating 
rate is given by 



n{t) = C^s{t)h{r)f^eom{r), 



(22) 



2^ where the radially dependent functions are nor- 
malized according to 



dr 47rr2 h{r) /geom(r) = 1 . (23) 



The geometrical factor is given by 

exp {-(r/rdisr)} . . 



4 



and £es{t) denotes the total 'effervescent' lumi- 
nosity at time t. 

In our continuous feedback model we assume 
that the luminosity of the AGN is proportional to 
the mass accretion onto the core, i.e. 

Ces{t) = eMcoreC^ (25) 

where the core is made up of gas whose entropy 
index lies below a minimal value, (Ti„i„. 

The accretion rate is calculated by a linearly 
weighted least-square regression of the function 
-^core(^), where a similar procedure as decribed 
in Sec. 2.3 is used with a smoothing timescale of 
20 Myr. In order to improve numerical stability 
we introduce a delay time of 2 Myr, i.e. the lu- 
minosity is proportional to the accretion rate of 
2 Myr ago. 

3. Results 

3.1. Self-regulated luminosities 

As our canonical model we take a cluster with 
a dark matter mass of Myir = 4 x 10^^ Mq and 
a concentration parameter c = 4. These values 
are chosen to resemble those of the cooling flow 
cluster AbcU 2052 (Blanton ct al. 2003). Wc start 
all our simulations from the initial temperature 
profile given by Eq. (2) and leave conduction and 
heating switched off until the cluster has formed a 
cold core. Only then we begin with the simulation 
of feedback. 

As the ICM cools and the central cooling times 
decrease, the accretion onto the core increases and 
so does, according to our prescription for feedback, 
the luminosity of the AGN. We will start by first 
describing the simulation without thermal conduc- 
tion. 

Fig. 1 shows the self-regulated luminosity of the 
effervescent heating as a function of time for var- 
ious values of the efficiency, e. After an initial 
phase which only lasts for a few hundred Myrs the 
luminosity converges to a value of around L ~ 10^'^ 
erg s~^ and remains steady for the entire length of 
the simulation, which is longer than 3 Gyrs. Runs 



with different efficiencies show, that the efferves- 
cent heating is in fact the regulating parameter. 
Even for quite different efficiencies that span two 
orders of magnitude the luminosities end up at 
quite similar values. 

The heating stalls the accretion onto the core 
as is shown in Fig. 2 which shows the accretion 
rate as a function of time. One can see how the 
accretion rate adjusts itself to a value just below 1 
Mq/yv after a few hundred Myrs. For comparison 
we also show the corresponding scenario with nei- 
ther heating nor conduction by the thick double- 
dashed line that is labelled by "all off" . It shows 
how quickly the accretion rate diverges in the ab- 
sence of heating. 

Many AGN are believed to be recurrent with 

duty cycles between 10 and a few 100 Myrs. This 
suggests that a non-linear feedback process is at 
work, i.e. that the luminosity is not linearly re- 
lated to the accretion rate as presumed here. How- 
ever, nothing is known aboTit the parameters of 
such a non-linear feedback mechanism. One could, 
for instance, introduce upper and lower thresholds 
for M, such that the AGN switches on if the ac- 
cretion rate lies above the upper threshold and is 
switched off when it falls below the lower thresh- 
old. Depending on the values for these thresholds 
one can reproduce sensible duty cycles, but since 
one has introduced two more parameters, this is no 
great achievement. Therefore, we decided to focus 
on a simple prescription for continuous feedback 
which introduces only a single parameter. Most 
importantly, our results depend only very weakly 
on this parameter. Furthermore, even if in nature 
feedback operates in some non-linear fashion, our 
results are still valid in a time-averaged (i.e. aver- 
aged over a few duty cycles) sense since the recur- 
rence times are small compared to the time scales 
considered here. Thus, if the AGN is active for 
only part of the time, the mechanical luminosities 
of the central AGN can be higher than the ~ 10^'^ 
erg s~^ found here. 
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3.2. The effect of conduction 

Next, we studied the effects of thermal conduc- 
tion on our models. Fig. 3 is the corresponding 
plot to Fig. 1 with thermal conduction (where we 
assumed / = 0.3). The solid lines in Figs. 1 and 3 
both correspond to the same efficiency of 3 x 10~^ 
and it is evident that the presence of thermal con- 
duction leads to decrease of the accretion rates 
and, consequently, also of the effervescent lumi- 
nosities. Initially, i.e. after about 500 Myrs, the 
luminosities are very similar at ~ 10'''* erg s~^ but 
at later times thermal conduction suppresses the 
luminosities to less than 10^^ erg s~^. Higher 
efficiencies reduce the power of the AGN even fur- 
ther, while efficiencies of around 3 x 10~^ result in 
a very steady luminosity. In any case, for a wide 
range of efficiencies we obtain cluster models that 
remain stable over long times and require 'efferves- 
cent' luminosities of not more than ~ 10^* erg s^*. 

Fig. 4 shows the volume heating and cooling 
rates of our simulations of clusters with masses 
of IQi'' Mq, 4 X 10*4 Mq and lO^^ Mg with an 
efficiency parameter of e = 3 x 10~^ and an as- 
sumed conductivity of 30% of the Spitzer conduc- 
tivity (/ = 0.3). The figure shows the curves at a 
time of 1 Gyr after the start of the feedback. The 
solid line shows the radial cooling profile, while the 
double-dashed line shows the effervescent heating 
rate. It is evident from the left panel that for the 
cluster with a mass of 10^^ Mq the radiative cool- 
ing is almost exactly balanced by the effervescent 
heating while conduction has very little effect on 
the energy budget. Also note, that thermal con- 
duction can, both, heat and cool the cluster as is 
shown by different lines in Fig. 4. In the outer 
regions and at the edge of the core, conduction 
leads to cooling while regions between a few kpc 
to ~ 100 kpc are heated. This has also been in- 
ferred from observations of M87 (Ghizzardi et al. 
2004). In the very center the heating due to con- 
duction drops to small values since the tempera- 
ture is very low. However, the relative importance 
of conduction increases with mass. This is not un- 
expected since the thermal conductivity depends 



strongly on temperature and temperature scales 
with mass. For the cluster with M = 4 x 10^** Mq 
conduction is already more important than effer- 
vescent heating for most regions within the cluster. 
Especially conduction can heat the region around 
the core and thus effectively suppresses the accre- 
tion. The figure shows that conduction alone al- 
most compensates the radiative cooling, particu- 
larly in regions where the bulk of X-ray emission 
comes from. This is in agreement with previous 
findings e.g. by Narayan & Medvedev (2001) and 
Fabian et al. (2002). 

3.3. Comparisons with observations 

While we do not intend to undertake a thorough 
modelling of the observed features of cooling flow 
clusters, it is still essential to verify whether our 
model can reproduce the gross features of typical 
cooling flow clusters. 

Fig. 5 shows the temperature profiles for a 
model with M^ir = 4 x Mq and c = 4.0 at dif- 
ferent times. Mass and concentration parameter 
are chosen to lie close to the values of Abell 2052, 
which is a well known cooling fiow cluster. The 
initial temperature profile is fairly flat as shown 
by the top curve in Fig. 5, but slowly the central 
temperature drops and after about 9 Gyrs a cold 
core has formed. The symbols in Fig. 5 denote the 
observed temperatures for A 2052 (Blanton et al. 
2003). We note that the data are consistent with 
our model at late times, i.e. after the formation of 
the core. Especially the slopes within 50 kpc agree 
nicely, whereas in the outer regions the observa- 
tions show higher temperatures than our model. 
This may be an indication that in the periphery 
other heating processes arc at work. Given the 
simplicity of our model, it may be imprudent to 
overstretch its similarity with observations but in 
any case it is encouraging that the profiles that we 
obtain are not disparate from observations. Again, 
we wish to stress that the temperature profiles are 
quite insensitive to the efficiency of the feedback. 

Finally, it is interesting to compare our accre- 
tion rates to the accretion rates inferred from the 
total x-ray luminosity according to the classical 



6 



cooling flow model, given by 

^cool-flow « \j-7^Lx{< Tc), (26) 

where is the cooling radius at which tc ~ 
Hq^ and Lx is the bolometric x-ray luminosity. 
Fig. 6 shows the inferred mass deposition rates 
for some of our simulated clusters. Comparing 
-^cool-flow to Mcorc rcvcals that the accretion rates 
inferred from the X-ray emission are considerably 
higher than the rates derived from our single-phase 
model. In the presence of heating, the simple 
equation for Mcooi-flow overestimates the accre- 
tion onto the core at least by a factor of 10. The 
reason is, of course, that the radiative losses are 
replenished by bubble heating and thermal con- 
duction, and thus unaccounted for in Eq. (26). 
Again, the variation of Mcooi-flow with e is rela- 
tively weak. To put this variation into perspective 
wc show how much more strongly Afcooi-fiow' de- 
pends on the concentration parameter (top curve 
in Fig. 6) when all other parameters remain the 
same. 

3.4. Summary 

We found that a simple prescription for feed- 
back leads to a self-regulated steady-state model 
for the ICM: For a wide range of efficiencies the 
self- regulated luminosities lie around L ~ 10^'' 
erg/ s which corresponds to a mass accretion rate 
of M ^ lM0/yr. Thermal conduction decreases 
the luminosity in a self-regulated model and is the 
more important the more massive the cluster is. 

The resulting profiles for temperature are con- 
sistent with observations of cooling flow clusters. 
In the presence of heating, the classical cooling 
flow model grossly overestimates the mass deposi- 
tion rates. 
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Fig. 4. — The volume heating and coohng rates for three models of different masses (from left to right the 
mass is 10^** Mq, 4 x 10^^ Mq and 10^^ Mq\ c = 4 for all models). The units on the vertical axis are erg 
cm~"^s^^. All three figures correspond to a time of 1 Gyr after the AGN has been switched on. The feedback 
parameter is e = 3 x 10""^. 



Cosmological model 

matter density Qo 0.3 

baryon density Q.b 0.04 

Hubble constant h = Hq/HI^^ 0.7 

Helium mass fraction Ijie 0.24 

Dark matter halo 

dark matter virial mass * -Mvir 10^^ Mq 

concentration parameter * c 4 

Effervescent energy feedback 

efficiency * e 3 x 1G~^ 

initial bubble size c?bubbie 20 kpc 

disruption radius raisr 500 kpc 

Conduction 

Spitzer fraction * /spitz 0.3 



Table 1: Parameters used for our standard simu- 
lation. To study other cluster models parameters 
marked with * are varied. 
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Fig. 1. — Luminosity of the effervescent heating as 
a function of time. The lines correspond to models 
with different efficiencies, e, The oscillations that 
appear at late times for very low luminosities are 
numerical artefacts caused by the finite resolution 
of the mass shells. 
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Fig. 2. Mass accretion rate onto the core as a 
function of time. The lines correspond to models 
with different efficiencies, e. The oscillations that 
appear at late times for very low luminosities are 
numerical artc^facts caused by the finite resolution 
of the mass shells. The thick double-dashed line 
that is labelled by "all off" is the accretion rate 
in a model without heating. It shows how quickly 
the mass accretion rate diverges in the absence of 
heating. 
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Fig. 3. — Same as Fig. 1 but with thermal con- 
duction. For comparison, the thin line denotes 
the volume integrated radiative luminosity of the 
ICM as a function of time. 




0.001 



r [Mpc] 



Fig. 5. — Temperature profiles for a cluster model 
with Mvir = 4 X 10" Mq and c = 4.0. Mass 
and concentration are chosen close to the values 
of Abcll 2052. which is a well known cooling flow 
cluster. The different lines correspond to different 
points in time with a new line drawn every 2 Gyr 
up to 16 Gyr. Symbols indicate the observed tem- 
perature profile (Blanton et al. 2003). 
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Fig. 6. — Mass deposition rates inferred according 
to the classical cooling flow model, i.e. Eq.(26). 
The dashed lines correspond to our canonical 
model (Mvir = 4 x lO^^M© and c = 4) with differ- 
ent efficiencies and conductivities, while the dot- 
ted fine corresponds to a model with £ = 10""^ but 
a higher concentration parameter (c = 7). 
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